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Abstract 

Energy transfer in the photo synthetic Fenna- 
Matthews-Olson (FMO) complex of Green Sul- 
fur Bacteria is studied numerically taking all three 
subunits (monomers) of the FMO trimer and the 
recently found eighth bacteriochlorophyll (BChl) 
molecule into account. The coupling to the non- 
Markovian environment is treated with a master 
equation derived from non-Markovian quantum 
state diffusion. When the excited state dynam- 
ics is initialized at site eight, which is believed to 
play an important role in receiving excitation from 
the main light harvesting antenna, we see a slow 
exponential-like decay of the excitation. This is 
in contrast to the oscillations and a relatively fast 
transfer that usually occurs when initialization at 
sites 1 or 6 is considered. We show that different 
values of the electronic transition energies found 
in the literature can lead to large differences in 
the transfer dynamics and may cause additional 
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suppression or enhancement of oscillations. 
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Figure 1: The Fenna-Matthews-Olson photosyn- 
thetic complex trimer, consisting of symmetry- 
equivalent monomers A, B, and C. Top: side view 
and labeling of the BChl molecules within one 
monomer. Bottom: top view. The figure was cre- 
ated using PyMOL^ and is based on the PDB entry 
3ENI. 

The Fenna-Matthews-Olson (FMO) complex 
bridges electronic excitation energy transfer from 
the chlorosome (the main light harvesting antenna 
complex) to the photo synthetic reaction center in 
green sulfur bacteria.-i^ It consists of three identi- 
cal subunits, which we denote as monomers A, B, 



and C (see [figure][l][]l ). Each of the monomers 



contains eight BChl molecules. Electronic ex- 
citation can be transferred between the BChl 
molecules via transition dipole-dipole interaction. 
Due to the availability of high resolution crystal 
structure a good overall understanding of the tran- 
sition energies of the chlorophyll molecules and 
their dipole-dipole couplings has been obtained, 
although there is still a large variance between 
values found in the literature (see Ref.- and in 
particular Ref.~). 

Regarding the arrangement of the FMO complex 
w.r.t. the chlorosomes, it is assumed that BChls 1, 



6, and 8 are located at the baseplate (which con- 
nects the chlorosomes and the FMO complex) 
and receive electronic excitation captured by the 
chlorosomes. In Ref. ^ evidence was found that 
the newly discovered eighth BChl directly obtains 
a significant part of the excitation handed over to 
the FMO complex. BChls 3 and 4 are located in 
the vicinity of the reaction center. This arrange- 
ment of the BChls is in accordance with the or- 
dering of their electronic transition energies. BChl 
8, BChl 1 and 6 which are located near the base- 
plate have large transition energies. Thus the FMO 
complex acts like a funnel that transfers energy 
from one side to the other, where BChls 3 and 
4 have low transition energies. Since the eighth 
BChl molecule was discovered only recently, it has 
not been considered in most of the previous studies 
(e.g. Refs.-i^"— ). Only in some recent studies the 
energy transfer has been treated taking the eighth 
BChl molecule into account.-^ 

In the present work we investigate theoretically 
excitation transfer in the full FMO trimer, focus- 
ing in particular on the role of BChl 8 and the de- 
pendence of the transfer on different sets of tran- 
sition energies. For our calculations we use a 
master equation which is derived from the non- 
Markovian quantum state diffusion (NMQSD) 
equation,—"— within the zeroth order functional 
expansion (ZOFE) approximation, ^^i^^i^^ jj^ ^ j.g_ 
cent publication we have shown (for the case of 
a single monomer subunit, taking only seven BChl 
molecules into account) that with this approach we 
obtain very good agreement with numerically ex- 
act results. This method, described in Ref.,^^ al- 
lows one to calculate the energy transfer in the full 
FMO trimer consisting of 24 coupled BChls within 
a few hours on a standard PC, taking coupling 
of the electronic excitation to the non-Markovian 
environment into account. Our master equation 
is derived from a stochastic Schrodinger equation 
by averaging over many trajectories.—'^ For in- 
dividual trajectories the coherence time might be 
substantially longer than for the average. Similar 
conclusions have been drawn recently. However, 
note, that the interpretation of such individual tra- 
jectories is far from trivial (see e.g.— >2^). 

Although the method enables us to treat compli- 
cated structured spectral densities for the exciton- 
vibrational coupling, in the present paper we 
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use a relatively simple spectral density shown in 
[figure] [2] []2[ which nevertheless leads basically 
to the same non-Markovian dynamics as the spec- 
tral density used in Ref.,-^^ as was shown in Ref.^*^ 
Important parameters are the interactions be- 
tween the chromophores and the site energies. For 
the interaction matrix elements between the BChl 
molecules we take values from Ref.,^ which are 



given in [table] [1][]1 for the interactions within 



one monomer and in [[table] [2][]2| for the interac- 
tions between the BChl molecules located in dif- 
ferent monomers. Note that these values and the 
values of the site energies used in the following 
are only according to a working model developed 
using biochemical and spectroscopic evidence. 

Following Ref.,- we assign BChl 8 to the 
monomers such that the strongly interacting neigh- 
bors BChl 1 and BChl 8 belong to the same 



monomer (see [figure][l][]l I and the coupling be- 
tween different monomers is relatively weak (< 
8 cm"^ see |[table][2][]2| ). Note that this assign- 
ment is different from the one of Ref.,^ where 
BChl 8 is assigned to the monomers based on the 
protein structure. 
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Figure 2: The spectral density employed, taken 
from Ref .^^ The reorganization energy of the spec- 
tral density is 35 cm^^ which is in the order of 
the interactions between the BChls and leads to 
a decoherence time of roughly 0.15 ps. The en- 
vironment correlation function corresponding to 
this spectral density is shown as inset (dashed line: 
imaginary part, solid line: real part). The correla- 
tion time is much shorter than the time range on 
which we consider the energy transfer. 

For the site energies of the BChls we con- 
sider values taken from Olbrich et. al. (OLB) in 



Ref.- and also values of Schmidt am Busch et. 
al. (SAB) in Ref., ^ which are both shown here 
in I [table] [3 ][] 3 1 Note the large differences be- 
tween the two sets of site energies. These two 
sets of transition energies for the BChls stem from 
recent publications and are based on the crystal 
structure. However, quite different strategies have 
been employed to obtain them. In Ref.- (SAB) 
quantum chemistry on BChl ground and excited 
states is combined with electrostatic calculations, 
where the whole protein is included in a very de- 
tailed fashion. In contrast, in Ref.^ (OLB) a time- 
dependent approach has been chosen in which 
many molecular dynamics ground state trajecto- 
ries have been run, calculating the electronic gap 
energy along the trajectory. This does not lead 
to a single transition energy but to a distribution 
for each BChl. This distribution was found to be 
asymmetric with tails towards higher energies and 
different for each BChl. In the present work we 
have taken the average positions of these distribu- 
tions to represent the electronic transition energy. 

Although there exist many sets of empirical pa- 
rameters in the literature (for a discussion see 
Ref.-i), we have chosen these two sets of values, 
since they are both obtained from a model taking 
the eighth BChl molecule into account. 

The calculated energy transfer in the FMO for 
the site energies of Ref.- is shown in [figure] [3] []3| 



and for those of Ref.- in [[figure] [4][]4[ The re- 
sults presented here are for a low temperatureQ of 
77 K for three different initial conditions. In the 
first row all excitation is located initially on BChl 1 
of monomer A. In the second row it is localized on 
BChl 6 and in the third row on BChl 80- 

The first thing to note is that for both sets of 
site energies and all three initial conditions there is 
only little transfer away from monomer A, which 
is due to the very small inter-monomer interactions 



(see [[table] [2] []2[ ). We have also investigated ini- 
tial conditions for which the excitation was delo- 
calized on the three monomers and found, as ex- 
pected, essentially the same behavior. 



'The values in OLB are based on a room temperature 
MD. However, we do not expect a strong temperature de- 
pendence for the average values. 

^Due to the rotational symmetry of the FMO trimer, this 
is equivalent to the case where the excitation would reside 
on monomer B or C. 
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Figure 3: Calculated transfer through the FMO 
trimer at T = 77 K for the energies of Ref. ^ (OLB). 
From left to right: Population on the BChls of 
monomer subunits A, B, C. First row: initial ex- 
citation on BChl Al. Second row: on BChl A6. 
Third row: on BChl A8. Note that the populations 
on monomer B and C are on a different scale than 
those for monomer A, because of the very small 
population that is transferred. 

We now take a closer look at the different ini- 
tial conditions. We first treat the transport along 
the FMO for the site energies of OLB and contrast 
them later to the transfer dynamics obtained for the 
site energies of SAB. 

One sees that for the initial excitation on BChl 1 
there are strong oscillations of the excitation be- 
tween BChls 1 and 2. This is due to the strong cou- 



Figure 4: As [[figure] [3] []3| but with the SAB en- 
ergies taken from Ref.- (see [[table][3][]3| ). 



pling between these two BChls (see [[table] [1 ][]![ ). 
Similar dynamics, but with much weaker oscilla- 
tions occur when the excitation is initialized on 



BChl 6 (second row of [figure] [3] [] 3 1. 
Different dynamics is found for the site energies 



of SAB shown in [figure] [4] []4 When starting on 



BChl 1 the oscillations between BChl 1 and 2 are 
even more pronounced and last longer than for the 
parameters of OLB. In the second row (starting on 
BChl 6) a stronger oscillation between BChl 6 and 
5 is now visible. A large difference in the transport 
between the two sets of site energies becomes ap- 
parent when considering the time dependence of 
the populations of BChl 3, which has the lowest 
site energy and thus in thermal equilibrium is ex- 
pected to have the highest population. While for 
the OLB energies there is already a saturation after 
about 0.6 ps , for the SAB site energies the popu- 



lations on BChl 3 are still growing at 1.0 ps. This 
strong difference can be attributed to the much 
lower energy of BChl 3 w.r.t. the energies of the 
other BChls for the SAB site energies compared 
to the OLB site energies. Note further that for the 
OLB site energies BChl 7 is detuned strongly from 
all other site energies and is thus excluded from 
the transfer. The effective spread of the site ener- 
gies, without BChl 7, is then only ~ 220 cm^^ and 
thus much smaller than in the SAB case, where it 
is ~ 500 cm^^ From this, one expects that the 
transfer is faster in the OLB case, in agreement 
with our observation. 

Remarkably, a completely different dynamics 
takes place when the eighth BChl is excited ini- 
tially (third row): In this case the excitation on 
BChl 8 decays very slowly with an exponential- 
like curve and no oscillations are present. This 
suppression of oscillations (see also Refs^>^) is 
due to the fact that the only relevant coupling of 
BChl 8 is to BChl 1 and that the population that 
is handed over to BChl 1 is immediately trans- 
ferred to BChl 2 due to the large coupling between 
BChl 1 and 2. Thus, no population can be accumu- 
lated on BChl 1 to initialize oscillations between 
BChl 1 and 2. A simple model is given in Ref.,— 
using the SAB site energies. In Ref.^^ the slow 
exponential-like decay is explained by the fact that 
the coupling of BChl 8 to BChl 1 is small com- 
pared to the energy gap between BChl 8 and 1 and 
the coupling to the environment. Indeed, for the 
SAB site energies the energy gap is 200 cm^^ and 
thus much larger than the coupling of 20 cm~^ 
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between BChl 8 and 1 . However, for the OLB en- 
ergies the gap is only 30 cm^^ i.e. comparable 
to the coupling between BChl 8 and 1 . We have 
found that for short times (< 0.1 ps) the dynam- 
ics is dominated by the (slow) coherent oscillation 
between BChl 8 and the strongly coupled dimer 
consisting of BChl 1 and 2. For later times deco- 
herence leads to the observed monotonic decay of 
the population. 

We also considered the transfer at a higher tem- 
perature of 300 K and found that it is similar to 
that at 77 K shown here, but the oscillations are 
less pronounced. When BChl 8 is excited initially, 
the transfer at 77 K and 300 K looks nearly the 
same. 

Furthermore, we also calculated the transfer for 
the energies and couplings of Ref.,^ obtained from 
the crystal structure instead of the MD simula- 
tions, and found almost identical transfer dynam- 
ics. 

Regarding the effects of the variations of the site 
energies, one should keep in mind that the val- 
ues of OLB stem from an atomistic model based 
on molecular dynamics simulations, which yield 
distributions for the site energies,-"^ which in the 
case of OLB show a pronounced asymmetry. In 
the present work we have used the mean values of 
these distributions as the site energies. We have 
also performed calculations using the maxima of 
the distributions and found qualitatively the same 
transfer dynamics as for the mean values. 

Beside the electronic transition energies and 
couplings also the interaction with vibrational 
modes plays an important role. While in our previ- 
ous work we have investigated structured spectral 
densities,— i^i^ here we have taken an often used 
featureless model spectral density. From atomistic 
simulations it is also possible to deduce spectral 
densities.—"^ While the one obtained in^^ has a 
much smaller effective reorganization energy (a 
measure for the strength of the coupling to the en- 
vironmentj^ than the spectral density used in the 
present work, the spectral density of OLB— has an 



^In Ref.— we defined the effective reorganization energy 
of the environment spectral density to quantify the part of the 
overall reorganization energy that is in the relevant frequency 
region, where the purely electronic dynamics couple to the 



effective reorganization energy which is an order 
of magnitude larger. We have performed calcula- 
tions for different effective reorganization energies 
by scaling the whole spectral density with a certain 
factor. As expected, we found that with decreas- 
ing/increasing reorganization energy the observed 
beatings become stronger/weaker. For an increase 



environment. For the spectral density of [figure] [2] []2 the 
effective reorganization energy is roughly 30 cm^^ 



of the original spectral density of [figure] [2] [] 2 by 
a factor of four (then the increased spectral density 
has roughly the same effective reorganization en- 
ergy as the averaged spectral density for BChl 1- 
6 obtained in Ref.^^), we observed that the beat- 
ings have vanished almost completely. Remark- 
ably, for this increased spectral density we found 
the population on BChl 3 after 1 ps to be about 
three times larger than the population for the orig- 
inal spectral density. This shows the strong in- 
fluence of the coupling to the environment on the 
transfer dynamics of the FMO complex. We have 
also performed calculations for a Markovian envi- 
ronment. Compared to the spectral density used in 
the present work, we found in the Markovian case 
that the population that arrives on BChl 3 after 1 ps 
is at least by roughly a factor of two lower - re- 
gardless of the coupling strength to the Markovian 
environment. 

In the present paper we considered the energy 
transfer through the full FMO trimer, where we 
also included the eighth BChl molecule. The 
values for the energies and couplings between 
the BChl molecules were taken from recent MD 
simulations and from the crystal structure (from 
Refs.^'^). For the calculation we used an efficient 
non-Markovian master equation (derived from the 
NMQSD). 

It was found, as expected, that due to the small 
couplings between the BChl molecules of different 
monomers (compared to the chromophore-protein 
coupling), there is only little exchange of elec- 
tronic excitation between the monomer units (this 
weak coupling limit has been discussed e.g. in 
Ref.-). This means that the three monomers act 
as independent transfer channels. Thus, the lim- 
itation to only a single isolated monomer unit of 
the FMO complex, as done in most of the previ- 
ous studies of the energy transfer in the FMO (e.g. 
Refs.—^~—), seems to be a reasonable approxima- 
tion considering the inter-monomer couplings and 
the short time scales adopted in the present work. 
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In the literature there exist many different sets 
of site energies of the BChls in the FMO (see 
e.g. Fig. 3 of Ref.^). Since these sets of ener- 
gies often differ quite strongly from each other, 
we have investigated the dependence of the trans- 
fer dynamics for two particular sets of energies 
(taken from Ref.- and Ref.-) and found large dif- 
ferences. In particular the amount of excitation ar- 
riving on BChl 3 (near the reaction center) may 
vary by more than 100% depending on which data 
set was taken. This shows that for further inves- 
tigations and interpretations of the transport (e.g. 
regarding an optimal efficiency) an accurate deter- 
mination of the energy landscape within the FMO 
complex is essential. 

Another crucial factor is the initial state, i.e. 
which BChl molecules are excited in the begin- 
ning. We considered different initial conditions, 
namely when the excitation is localized either on 
BChl 1, 6, or 8. This choice of the initial state is 
motivated by the assumption that BChl 1, 6, or 8 
have the highest probability to obtain the excita- 
tion from the chlorosome antenna structure. Re- 
markably, when initially only BChl 8 is excited, 
it was found that a slow exponential-like transfer 
away from BChl 8 takes place. This is in con- 
trast to the oscillations found when starting on 
BChl 1 and the faster transfer away from BChl 1 
or 6. These oscillations have been the subject of 
study of many research groups, both theoretically 
and experimentally. Ultrafast experiments where 
BChl 8 is present as well might reveal this strong 
difference in the energy transfer dynamics away 
from it. Furthermore, we found that the influence 
of variations of different parameters on the trans- 
fer strongly depends on the initial state. Here, 
the parameters we varied include electronic ener- 
gies, coupling energies between BChl molecules, 
or the coupling to the environment. Because of 
this strong dependence on the initial condition, in 
further investigations it would be important to get 
more detailed information about the arrangement 
of the FMO w.r.t. the baseplate and the chloro- 
somes and learn more about the transfer of exci- 
tation from the chlorosomes to the FMO complex. 

Let us briefly compare our findings with those 
of Refs.,->^ from which the couplings and ener- 
gies of the BChls were taken. In Ref.^ Olbrich 
et. al. found an even slower decay of the excita- 



tion on BChl 8 than in the present work. As in 
the present paper, they also observed a much faster 
transport away from BChl 1 or 6 compared to 
BChl 8, but no oscillations, which can be attributed 
to the much higher (effective) reorganization en- 
ergy of their spectral density. Schmidt am Busch 
et. al. found a considerably faster decay of the ex- 
citation on BChl 8. ^ This can be explained by the 
fact that in their calculation the strong coupling of 
BChl 8 to BChl 1 is almost twice as large as the 
respective coupling in Ref.,- which was taken in 
the present work. Note that in the calculations of 
Ref. ^ the initial state is a superposition of excita- 
tion on BChls 8, 1, and 2. 

We have also investigated various initial states 
where the excitation is coherently distributed over 
different BChls (which may be located on different 
monomers) to investigate the possibility of phase 
directed transport. For all cases considered the 
initial phase had a marginal effect on the transport 
dynamics. However the initial population strongly 
influences the transport dynamics. Thus in further 
studies to avoid the arbitrary choice of the initial 
state and to take the flow off of excitation into ac- 
count, also the chlorosome (at least the part close 
to the FMO) and the reaction center (which plays 
a significant role) should be included into a holis- 
tic simulation. Efforts towards this goal are being 
made in our research groups. 
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Table 1 : Intra-monomer couplings from Ref for 
the trimer structure in cm ^ obtained using MD 
simulations. Values greater than 10 cm^^ are high- 
lighted in bold face. 
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-8.7 



Table 2: Inter-monomer couplings from Ref.- for 
the trimer structure in cm ^ obtained using MD 
simulations. Upper triangle: couplings between 
monomer units A-B, B-C and C-A. Lower trian- 
gle: couplings between monomer units A-C, B-A 
and C-B. 
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Table 3: Site energies of BChl 1-8 on each 
monomer. The values are given in cm^^ The 
respective lowest value is set to zero (for the cal- 
culation of the transfer only the energy differences 
are relevant). Note that for the values of Olbrich et 
al. the average positions from Table 1 of Ref.^ are 
taken. 
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